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J ' Hils  report  proposes  robust  detection  theoretic  approaches  to  HF  signal  reception.  Most  of 
these  approaches  require  modifying  recent  robust  detection  theory  results  to  the  specifics  of  the 
HF  communication  problem.  One  approach  recommended  la  direct  and  avoids  this  modification 
requirement;  however,  it  demands  more  original  effort.  Hie  report  contains  a survey  of  recent 
robust  detection  theory  research  and  discusses  the  applicability  of  the  research  results  to  HF 
additive  interference  models;  It  concentrates  on  the  known-signal  case  to  highlight  the  additive 
interference  problem. 
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ON  APPROACHES  TO  ROBUST  DETECTION  FOR  HF  COMMUNICATIONS 


INTRODUCTION 

Simply  stated,  robust  detection  capitalizes  on  all  the,  a priori,  noise 
model  information  available  to  the  designer  and  guarantees  a detector  performance 
level  against  all  noise  models  described  by  that  Information.  It  Is  suboptlmal 
for  each  noise  model  allowed  by  the  a priori  Information,  but  generally 
Insensitive  to  changes  from  model  to  model.  Suppose  you  have  substantial,  but 
Incomplete,  knowledge  of  the  additive  noise  process.  Optimal  detection  would 
be  Impossible  since  It  requires  complete  specification  of  the  noise  probability 
density-— you  probably  would  design  for  a noise  that  occurs  occasionally  or 
not  at  all.  Nonparametrlc  detection  would  be  too  conservative  since  It  assumes 
Insufficient  knowledge  for  limiting  the  noise  probability  density  functions  to 
a class  characterizable  by  a finite  number  of  parameters  — you  would  protect 
against  noise  that  you  know  will  not  occur.  For  example,  a robust  detector  for 
nearly-Gausslan  noise  ([16])  would  out-perform  a nonparametrlc  detector  (e.g., 
a sign  detector)  In  that  noise  environment  and  would  be  slightly  out-performed 
by  the  optimal  detector  In  Gaussian  noise.  Conceptually,  robust  detection  Is 
a game-theoretic  ([10],  [20]),  yet  pragmatic,  approach  to  detecting  and 
distinguishing  signals  In  an  Incompletely- known  noise  environment. 

Consider  the  robust  detection  problem  as  a minimax  or  game  theory  problem: 
Interpret  the  robust  (minimax)  detector  as  the  best  detector  strategy  (player  1) 
against  any  of  nature's  (player  2)  possible  noise  model  choices  within  the  a 
priori  class  of  choices.  More  abstractly,  consider  the  following: 

Note.  Manuscript  submitted  November  2,  1978. 
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PROBLEM:  Determine  the  detector  $°  and  least-favorable  Interference 

probability  density  vector  7°  = (fjj,  f? f°_j)  for  the  M-ary 

hypothesis  testing  problem, 

Hr*  s1  (t)  + WjU)  1 = 0 M-l,  0<t<T, 

such  that 

C(*°,  7°)  = min  max  C($,  7) 

HO  f«F 

where  C($,  7)  Is  the  cost  Incurred  by  the  detector  $ when  the 
Interference  Is  defined  by  7,  D Is  the  class  of  admissible  detectors, 
and  F is  the  class  of  Interference  probability  density  vectors. 

The  robust  (minimax)  detector  can  be  Interpreted  as  the  best  (with  respect 
to  C)  choice  from  0 for  a class  of  possible  descriptions  F.  Thus,  the  optimal 
(til],  [21],  [24],  [28])  and  the  nonparametrlc  ([4],  [27])  detection  problems 
are  the  opposite  extremes  (F  contains  one  vector  and  F contains  all  possible 
density  vectors,  respectively)  of  the  robust  detection  problem  and  can,  in  fact, 
be  encompassed  by  It. 

This  report  proposes  robust  detection  approaches  to  HF  signal  reception. 

Most  of  these  approaches  require  modifying  recent  robust  detection  theory 
results  to  the  specifics  of  the  HF  conmunlcatlon  problem.  One  approach 
recommended  Is  direct  and  avoids  this  modification;  however.  It  demands  more 
original  effort.  The  report  contains  a survey  of  recent  robust  detection  theory 
research  and  discusses  the  applicability  of  the  research  results  to  HF  additive 
Interference  sources  and  models;  It  concentrates  on  the  known-signal  case  to 
highlight  the  additive  Interference  problem.  Remember,  the  unknown-signal  case, 
which  would  be  appropriate  for  treating  fading  and  multipath  effects.  Is 
essentially  the  known-signal  case  either  appropriately  conditioned  (probabilistically) 
on  a random  signal  variable  or  assigned  to  some  least  favorable  known-signal  value. 

A subsequent  report  will  treat  the  unknown-signal  case  using  this  probabilistic 
conditioning  technique. 
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RECENT  RESULTS  IN  ROBUST  DETECTION  THEORY 


General  Contaminated  Noise  Model  — Hypothesis  Dependent 

The  current  Interest  in  robust  detection  began  with  the  work  by 
Huber  ([13])  on  robust  probability  ratio  tests.  He  was  searching  for  ratio 
tests  that  were  insensitive  to  a few  bad  observations  (extremity  data  points): 
so  that  those  few  bad  observations  could  not  override  the  majority  decision. 

A general  binary-hypothesis  independent-measurements  model  was  assumed  with 
different  contaminated-nominal  sample  distributions  for  each  hypothesis; 
in  other  words,  the  actual  distributions  were  known  up  to  some  neighborhood 
of  the  nominals.  The  class  of  possible  distribution  measures  for  each 
hypothesis  was  given  by 

Pi  = (Q1  |Qi  ■ (1  - ei)P1  + e^,  H^H)  i = 0,  1 
where  0 < ci  < 1 are  fixed  numbers,  and  H denotes  the  class  of  all  probability 
measures  on  some  given  measurable  space.  This  formulation  allowed  hypothesis- 
dependent  contamination  and,  consequently,  permits  interference  generated 
by  an  intelligent  adversary. 

Huber  obtained  the  worst-case  distribution  functions  and  robust 
(minimax)  ratio  tests  with  respect  to  three  standard  performance  measures. 

The  performance  measures  are  based  on  a risk  function  ( R(Q-j . 4>)  and  given  as 

i)  max  R(Q1t  $) 

ii)  R(Qj,  <|>)  subject  to  R(Qq,  $)  <_  a (Neyman-Pearson) 

111)  X0R(q0,  <J>)  + XjRfQj,  4>)  (Bayes) 

His  results,  which  are  valid  for  both  fixed  sample  size  and  sequential 
problems,  require  small  enough  to  prevent  overlapping  of  the  distribution 
classes  (preventing  one  distribution  for  both  hypotheses).  The  worst-case 
(least-favorable)  distribution  pairs  (Q0,  Qj)  generate  probability  ratios 
that  are  censored  versions  of  the  nominal  probability  ratios: 


■%- .. . 


be 


Pj  (x)/p0(x)  < c ' 
q1(x)/qQ(x)  = /bp^xj/p^x)  c'<  P1(x)/pQ(x)  < c" 

(be"  P^xJ/PqU)  > c" 

where  0 < c'  < c " < °°  are  chosen  so  that  qj  and  q0  are  probability  densities, 
and  b = (l-ejJ/d-eg) . Robust  ratio  tests  <f>°  are  acquired  via  the  Neyman- 
Pearson  Lemma,  and  the  above  probability  ratios.  Thus  the  robust 
ratio  tests  for  these  performance  measures  are  censored  versions  of  the 
optimal  ratio  tests  for  the  nominal  densities  — an  asymmetrical  soft- 
limiter  is  placed  within  the  optimal  detector  for  the  nominal  densities. 

Huber  also  obtains  equivalent  results  for  the  class  of  possible  distribution 
measures  described  by 

Pi  = tQi  I llQi  - Pill  ±e}  i = 0.  1 

where  ||  • ||  denotes  total  variation. 

Contaminated  Gaussian  Noise  Model  - Hypothesis  Independent 

Martin  and  Schwartz  ([16])  modified  Huber's  results  and  applied 
them  to  the  contaminated- normal  distribution  model  for  the  fixed-sample  size, 
signal  detection  problem.  They  sought  robust  detectors  of  known 
signals  (possibly  time-varying)  in  nearly-Gausslan  1.1. d.  noise.  The 

P 

possible  noise  distributions  were  assumed  to  be  of  the  mixture  model  form 
F(x)  = (l-e)'t(x)  + eH(x)  for  0 £ e < 1, 

where  4>  is  the  unit  normal  distribution,  H is  an  arbitrary  distribution,  and 
e is  small.  Both  regular  and  small-signal  cases  were  treated.  Some  atten- 
tion was  given  to  the  Incoherent  signal  detection  problem. 

For  the  nearly-Gausslan  noise  model  x^  * 0Sf  + n^ , 1 = 1,  2,  . . , N, 
Martin  and  Schwartz  found  a time-varying  correlator-limiter  as  the 
robust  (minimax)  solution  for  given  e and  0 > 0c/min^|s1|  with  respect  to 
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the  Neyman-Pearson  performance  measure.  Explicit  expressions  for  the  worst- 
case  (least  favorable)  densities  were  obtained.  Table  1 displays  representative 
values  of  e and  e^.  Note,  the  minimum  received  signal-to-noise  ratio 
(SNR)  can  be  derived  from  0e  via  the  equation  (SNR)min  = 20  log^^!p£-j.  For 
a constant  signal  (s^  = c),  the  minimax  solution  is  a 1 imi ter -detector  since 
correlation  is  unnecessary. 

For  the  small-signal  problem  with  nearly-Gaussian  white  noise,  the 
limiter-correlator  was  found  the  asymptotical ly-robust  solution  to 
the  local  ( small -signal ) Neyman-Pearson  problem.  Explicit  expressions  for 
the  worst-case  symmetric  densities  were  obtained.  This  result  is  valid  for 
false-alarm  probability  a ^a(e),  depending  on  given  e,  and  for  symmetric  contamina- 
tion densities  satisfying  a regularity  condition.  It  is  valid  also  for  a 
sample  size  N = 1 and  all  0 £ a £0.5.  The  authors  conjectured  the  validity 
of  this  result  for  a restricted  range  of  a,  which  depends  on  N.  Table  2 
gives  typical  values  of  e and  a(e).  The  limiter  in  this  case  is  a 
symnetrical  soft-limiter  with  break  points  K (see  table  2).  Again  with 
constant  signal,  the  resulting  robust  detector  is  the  limiter  detector. 

Martin  and  Schwartz  applied  the  soft  limiter  to  the  envelope  sum 
detector  which  is  the  optimal  small -signal  detector  for  an  Incoherent  detection 
of  a pulse  train.  The  resulting  limiter-envelope  sum  detector  was  shown  to 
exhibit  (for  the  examples  chosen)  the  same  degree  of  robustness  In  terms  of 
asymptotic  relative  efficiency  (ARE)  as  does  the  limiter-correlator. 


Table 

1 

Table  2 

*[!sec]  ■ >5(1 

■e)"1 

e 

9e 

20  log  6e 

c 

a(e) 

K 

0.01 

0.025 

-32.04  (dB) 

0.01 

0.163 

1.95 

0.012 

0.03 

-30.46 

0.02 

0.166 

1.72 

0.02 

0.05 

-26.02 

0.05 

0.174 

1.40 

0.04 

0.10 

-20.0 

0.10 

0.184 

1.14 

0.055 

0.15 

-16.49 

0.06 

0.16 

-15.92 

0.10 

0.28 

-11.06 

0.138 

0.40 

-7.96 

Contaminated  Non-Gaussian  Noise  Model  — Hypothesis  Independent 

Kassam  and  Thomas  ([14])  generalized  the  results  of  Martin  and  Schwartz 
for  asymptotically-robust  local  detection  of  a known  time-varying  signal  by 
using  a contaminated  1.1. d.  noise  model  with  a non-Gaussian  nominal.  The 
results  are  shown  valid  for:  nominal  density  functions  that  are  symmetric, 
strongly  unlmodal  on  Its  support,  and  have  an  absolutely  continuous  derivative 
within  Its  support;  contamination  density  functions  that  are  bounded  and 
symmetric;  and  false  alarm  probabilities  a > a(rm),  where  rm  Is  a parameter 
depending  on  c.  Kassam  and  Thomas  used  the  methods  of  Martin  and  Schwartz 
and  Huber  to  obtain  the  worst-case  noise  density  and  the  asymptotically- 
robust  local  detector  that  Is  the  optimal  local  detector  for  that  density. 

The  parameter  rm  depends  also  on  the  nonlinearity  function  (ZNL)  of  the 
optimal  local  detector. 

The  results  were  applied  explicitly  to  zero-mean,  general  1 zed-Gauss Ian 
nominal  densities  that  are  parameterized  by  their  rates  of  exponential  decay. 
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For  this  class  of  nominal  densities,  a ( r m ) , for  a Gaussian  nominal,  is  a 
lower  bound  on  a than  the  lower  bound  obtained  by  Martin  and  Schwartz  (see 
table  3).  The  general  asymptotically- robust  local  detector  solution  for 
this  class  of  nominal  densities  corroborates  the  limiter-correlator 
detector  for  contaminated  Gaussian  noise  and  Implies  the  robustness  of 
the  sign  detector  for  contaminated  double-exponential  noise.  For  weaker 
assumptions  on  the  contamination  densities  (zero  median  and  continuous 
at  the  origin),  the  authors  show  the  sign  detector  to  be  the  asymptotically 
robust  local  detector  for  all  0 <_  a < 1. 

Table  3 

(Gaussian  nominal ) 


e 

«(rm) 

0.02 

0.103 

0.04 

0.093 

0.06 

0.086 

0.08 

0.081 

0.10 

0.076 

Non-Contaminated-Noise  Model  — Hypothesis  Independent 

El-Sawy  and  VandeLinde  ([8],  [9])  used  a non-contamination  (mixture) 
model  for  the  noise:  to  eliminate  the  limitations  of  earlier  work  caused  by 
small  e.  They  used  also  a restricted  class  of  detectors  to  reduce  the  number 
of  detectors  considered  for  their  asymptotically- robust  detection  problem. 

The  binary  hypothesis  model  for  this  problem  takes  the  form: 

H0:  x1  * w1 

Hl:  X1  “ 9lsi  + , i * 1,  . . . , N 

where  {x^}  Is  the  sequence  of  observations,  6^  > 0,  (s-j } Is  a sequence  of 
known  constants  with  | s -f | < <*>, 


1 1m  N'1^  s2  = c, 

If*®  1=1  1 

and  {wi}  Is  a sequence  of  1.1. d.  random  variables  with  a syimietric  density  f. 

A generalized  version  of  Huber's  ([12])  M-estimate  of  location  was  introduced 

via  a class  of  functions  (*  satisfying  certain  smoothness  assumptions.  The 

generalized  M-estimates  0^  — the  admissible  test  statistics  — are  defined 

as  the  values  of  e which  minimize 
N 

L L(x1  - es^,  Le<J 
1*1 

The  restricted  class  of  detectors  deO  considered  are  just  the  threshold 

tests  obtained  from  the  generalized  M-estimate  6n  test  statistics.  The 

authors  solved  a constrained  asymptotic  (in  N)  maximln  problem  with  respect 

to  the  power  function  6^(0, |f).  They  derived  the  saddle-point  pair 

(d^*.  f*).  where  f*  is  the  minimum  Fisher-information  symmetric 

density  and  d^  is  a threshold-test  detector  for  the  generalized 

* * 

M-estimate  test  statistic  derived  from  L = -log  f . The  threshold  is 
chosen  to  yield  the  desired  false  alarm  rate  a = ^(Olf).  Equivalent 
results  are  obtained  for  the  local  detection  problem 
(subject  to  a lowerbound  on  a)  and  for  the  detection  problem  (derived  with 
(si } constant)  with  a class  of  p-point  noise  distributions  defined  as 
{f :/_!  f(x)dx  = p,  f symmetric}. 

El-Sawy  and  VandeLinde  claimed  also  that  the  results  may  be  used  also  to  solve 
other  hypothesis  testing  problems,  e.g.,  detection  problems  for  the  minimum 
probablllty-of- error  criterion. 

A Robust  Detection  Problem  for  HF  Comnuni cations 

The  robust  detection  problem  for  HF  conmunlcatlons  is  derived 
from  the  general  robust  detection  problem  by  being  more  specific  about  the 
cost  function  and  the  HF-nature  of  the  model  variables.  Using  a generalized 
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binary  hypothesis  model  with  the  probability- of- error  cost  function,  the 
robust  detection  problem  for  HF  communications  becomes  one  of  des.gning  a 
robust  test  for: 

: x(t)  = Si  (t ) + w-j  ( t ) , i = 0,  1;  0 <_  t ^ T 
where  s^(t),  with  | s -e ( t ) | < «>,  is  the  transmitted  waveform  after  the  channel 
filtering  effects  (front-end  filters,  channel  attenuation,  multipath,  etc.) 
and  w.j(t)  is  the  additive  interference  process,  with  | Wi  ( t ) j < °°,  generated 
by  an  additive  combination  of  atmospheric  noise,  thermal  noise,  man-made 
noise,  other-user  interference,  and  intelligent-adversary  interference.  In 
characterizing  the  additive  interference  process  w^(t)  as  an  element  from  a 
class  of  possible  stochastic  processes  at  the  minimum,  we  must  admit  the 
common  factors  of  the  accepted  (proposed)  models  for  the  different  types 
of  interference. 

Atmospheric  noise  has  been  modeled  in  several  ways  as  an  impulsive, 
non-stationary,  correlated  random  process  ([5]-[7],  [23],  [26])  usually  as 
a combination  of  a low-level,  high  probability  component  whose  envelope  is 
Rayleigh  distributed  and  of  a high-level,  low  probability  component  whose 
envelope  is  log-normal  distributed.  The  most  recent  models  ([23],  [26]) 
assume  the  impulsive  interference  waveforms  are  emitted  according  to  the 
Poisson  distribution  in  time  and  are  superimposed  on  an  additive,  independent 
Gaussian  background -noise  process.  These  impulsive  waveforms  may  be  modeled 
as  the  output  of  time-varying,  possibly  stochastically  described,  linear 
systems  excited  by  Poisson-distributed  Impulses.  A general  atmospheric 
noise  model  results  described  as 

wA(t)  = Wg(t)  + Wj(t) 

-9- 


where  Wg(t)  is  zero-mean,  independent  Gaussian  process,  and 

N(t) 

Wjlt)  = 22  u 
i=l 

with  hi  a possible  impulse  response,  ui  a finite,  amplitude  random  variable, 
t ^ the  point  process  variable,  and  N(t)  a counting  process  ([25]).  Thermal 
noise  that  is  white  Gaussian  and  man-made  noise  that  is  impulsive  are  of 
secondary  concern  for  HF  communications  ( 6);  nevertheless,  they  can  be 
adequately  handled  by  this  atmospheric  noise  model.  Hence,  no  further 
reference  will  be  made  to  these  two  sources  of  interference. 


Other-user  interference  is  attributed  to  randomly- occuring 
modulated  waveforms  of  various  durations  and  strengtns.  These  waveforms 
are  generally  confined  to  a 3 kHz  bandwidth  (’6.) — there  are  some  exceptions 
such  as  spread-spectrum  modulation  and  coding.  The  other-user  interference 
waveform  WQy(t)  also  can  be  modeled  as  a random  process,  possibly  in  the 
frequency  domain,  since  most  available  information  is  from  spectra 
measurements.  For  example,  consider  the  random  point-process  model  of 
figure  1. 


Figure  1 


This  is  a frequency  domain  description  of  other-user  interference  given  by 

M W f 

w0u(n  = L E U,  / Hk(f-X)6(,\-f1  )dX, 

0U  k=l  ik=l  k -oo  k 

which  implies, 

M Nk(f) 

wou(t)  - E E -i  khk(t)ej2  V . 

k=l  ik=l 

H, (f)  is  one  of  the  M possible  transmitted  transforms,  u.  is  the 

’k 

finite,  signal -ampl itude  random  variable  for  the  ikth  carrier,  fj  is  the 
ikth  carrier-frequency  point-process  variable,  and  Nk(f)  is  a counting  process. 

Intelligent-adversary  interference  (jamning)  can  take  any  form 
within  the  physical  constraints  Imposed  by  the  Interference  generating 
equipment.  These  constraints  usually  are  just  upper  bounds  on  peak  and 
average  power,  and  bounds  on  bandwidth.  This  characterization  admits 
stochastic  (e.g.,  wideband  Gaussian  noise)  or  purely  deterministic  (e.g.,  CW) 
interference  that  is  ideally  signal-waveform  dependent  for  maximum  effectiveness. 
For  signa 1 -waveform  dependence,  the  interferer  requires  knowledge  of  the  wave- 
form transmitted;  therefore,  we  assume  the  interferer  knows,  a priori,  the  set 
of  waveforms  transmitted,  but  not  the  particular  waveform  transmitted.  With 
these  assumptions,  the  intelligent-adversary  interference  can  be  modeled  as  a 
stochastic  process  WjA(t)  described  by  some  probability  density  function  (pdf), 
f(wIA(t),  t)ef,  where  fe  f implies: 


(a) 

f 

is  a first-order 

pdf 

(b) 

E 

8 

V 

o 

vl 

r-*-» 

-*-> 

<VI 

* 

0 < t < T 

(c) 

f( 

x ,t  ) « 0 

for  |x|  > c2  , 0 <_  t < T 
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APPLICABILITY  OF  RECENT  RESULTS 

The  recent  results  on  robust  detection  theory,  all  using  discrete-time 
models,  are  applicable  In  only  a limited  sense  to  the  robust  detection  problem 
for  HF  communications  as  discussed  earlier.  The  main  difficulty  Is  the  assumption 
of  1.1. d.  Interference  random  variables  which  Is  necessary  to  achieve  the 
results  primarily  based  on  asymptotic  normality  or  multiplicative  probability- 
ratio  test-statistics.  The  mixture  model  used  In  several  treatments  also 
presents  some  difficulty.  It  assumes  a time-sharing  of  statistical  descriptions: 
one  required  to  occur  only  a small  percentage  of  (e)  of  the  time  and  the  other 
a nominal  required  to  be  completely  known  statistically.  Finally,  the  require- 
ment of  symmetrical  Interference  densities  for  validity  of  some  results  Is 
possibly  a very  restrictive  constraint  on  interference  characterization  for 
this  problem. 

This  assumption  of  1.1. d.  Interference  random  variables  certainly  Is  not 
met  by  sampled  values  of  Interference  from  atmospheric  noise  or  other-user 
Interference;  and  It  may  be  a costly  limitation  for  Intelligent -adversary 
Interference  although  that  has  to  be  determined.  Although  techniques  for 
accommodating  the  non-l.l.d.  nature  (primarily  the  dependence)  of  the  HF 
Interference  are  desirable  th'»  Independence  assumption  may  be  used  to  determine 
bounds  on  detector  performance.  The  general  argument  Is  that  dependence 
(correlation)  among  samples  provides  additional  Information  for  designing 
counter-strategies,  thus,  the  assumption  of  Independence  leads  to  bounds  on 
performance  for  the  dependent  case.  In  addition,  the  Identically-distributed 
assumption  may  not  be  required  for  validity  results  previously  mentioned  since 
It  Is  not  necessary  for  asumptotlc  normality  ([1],  ch.  8)  nor  required  for 

• ] 

the  canonical  form  of  the  probability  ratio  tests  ((20)). 
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The  mixture  models  assumed  for  the  Interference  density  require  complete 
statistical  knowledge  of  the  nominal  density;  yet  for  the  HF  Interference 
discussed,  all  of  the  components  — atmospheric  noise,  other-user  Interference, 
and  intelligent-adversary  Interference  (when  applied) — are  not  known  exactly. 
Thus,  the  reported  results  for  this  mixture  model  could  not  be  applied  directly 
to  the  general  HF  coitmunlcatlon  problem  as  formulated  here.  You  could  Interpret 
the  mixture  model,  with  Its  time-sharing  Implication,  as:  (1)  a high  probability 
(1-e)  nominal  density  describing  all  but  the  Impulsive  noise,  and  a low 
probability  (e)  contaminated  density  describing  the  nominal  plus  the  impulsive 
noise  ([17]);  or,  (2)  a high  probability  nominal  Including  all  but  the  Intelligent- 
adversary  Interference  ([26])  and  a low  probability  contaminated  density 
describing  the  nominal  plus  Intelligent-adversary  Interference.  In  other  words, 
the  contaminating  Interference  component  and  the  other  noise  components  would 
be  assumed  to  occur  with  an  average  rate  of  e.  The  mixture-model  results  may 
be  applied  to  several  typical  nominal  components;  the  outcome  may  be  Instructive, 
typical  densities  for  atmospheric  noise,  other-user  Interference,  Intelligent- 
adversary  Interference,  or  any  combination  would  be  the  candidates  for  the 
nomlnals.  The  robust  detectors  obtained  from  the  different  nomlnals  could 
then  be  compared  for  similarities  that  could  be  exploited  later.  Another 
difficulty  In  applying  Huber's  mixture-model  results  to  the  HF  problem  Is  the 
assumption  of  a signal -dependent  Interference  capability  that  Is  not  assumed 
In  the  HF  communication  problem  formulated  here.  If  possible.  It  may  be 
Instructive  to  consider  detectors  for  signal-independent,  contaminating 
Interference  derived  from  a random  or  non-ramdom  convex  combination  of  the 
least-favorable  signal -dependent  densities. 

In  the  problem  formulated  here  for  HF  conmunlcatlons,  the  cost  function 
prescribed  Is  the  probability  of  error  (Bayes)  performance  measure.  However, 

-13- 


— - ^ 


only  Huber  ([13])  considered  this  cost  function  explicitly.  El-Sawy  and 
VandeLInde  Indicated  their  techniques  are  applicable  to  the  probability  of 
error  cost  function.  The  cost  function  considered  In  all  of  the  research 
discussed  Is  the  one  appropriate  to  radar  problems,  namely,  the  probability 
of  detection  failure  (miss  probability  1-e).  However,  since  the  detectors 
chosen  with  respect  to  this  radar  cost  function  are  generally  required  to 
satisfy  a probability  of  false-alarm  (a)  constraint  (the  Neyman-Pearson 
problem),  the  solution  for  the  communication  (Bayes)  problem  can  be  arrived 
at  by  solving  a series  of  radar  problems.  This  follows  from  the  fact  that 
the  probability  of  error  Is  a weighted  sum  of  the  miss  and  false-alarm 
probabilities  and,  consequently,  the  minimum  probablllty-of-error  detector  Is 
also  a minimum  probablllty-of-mlss  detector  achieving  a particular  a ([11], 
sec.  3.2).  Unfortunately,  the  results  reviewed  above  for  the  Interference- 
mixture-model  Neyman-Pearson  problem  are  valid  only  for  false-alarm  probabilities 
larger  than  a given  lower  bound  a(c)  and,  therefore,  may  not  be  applicable  to 
the  HF  communication  problem  formulated  above.  Unless  the  Bayes  detector  can 
be  found  among  those  Neyman-Pearson  detectors  satisfying  the  lower  bound 
constraint  a(c),  the  results  for  the  mixture  model  problems  cannot  be  applied. 

All  attention  must  then  be  g*ven  either  to  applying  the  results  ol  Huber  and 
of  El-Sawy  and  VandeLInde  or  to  solving  directly  the  robust  detection  problem 
for  HF  communications. 

The  requirement  of  symmetry  for  the  Interference  densities  may  be  a 
very  restrictive  limitation  for  this  HF  problem.  This  synmetry  of  densities 
Implies  zero-mean  Interference  and  certainly  eliminates  other-user  Inter- 
ference. Moreover,  It  Is  not  known  at  the  outset  If  restricting  Intelligent- 
adversary  Interference  to  being  zero-mean  Is  severly  limiting  or  not 
for  this  problem.  For  example;  Blachman  ([2],  [3]  ch.  9)  has  shown  that 
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zero-mean,  white  Gaussian  noise  Is  the  best  Intelligent-adversary  Interference 
for  minimizing  the  channel  capacity  seen  by  a comnunlcator;  and.  It  can  be 
shown  ([19],  [20])  that  deterministic  signals  are  the  best  Interference  for 
minimizing  the  correlation  In  correlator  receivers.  Thus,  It  appears  the 
symnetry  requirement  possibly  eliminates  or  sharply  curtails  the  two  main 
components  of  the  Interference  characterized  In  the  HF  problem  formulated  above. 
RECOMMENDATIONS  FOR  APPROACHES 

From  the  review  of  the  recent  results  on  robust  detection  and  the  discussion 
of  their  applicability  to  the  robust  detection  for  HF  communications  problem, 
recommendations  can  be  made  for  approaching  this  problem.*  The  first  approach 
recommended  Is  based  on  Huber's  ([13])  formulation  and  results.  For  the  second 
approach,  a direct  solution  attempt  Is  recommended.  In  both  cases,  the  concern 
Is  for  the  binary  detection  problem  where  both  signal  values  are  non-zero.  In 
addition,  both  approaches  are  based  on  a discrete-time  model  In  anticipation 
of  digital  processing. 

Approach  1 

In  the  first  approach;  the  mixture  model,  the  lack  of  symmetry  In  the 
Interference  densities,  and  the  treatment  of  the  finite  sample  problem  — instead 
of  the  asymptotic  problem — are  the  motivating  factors  for  utilizing  Huber's  results. 
The  mixture  model  for  the  Interference  densities  will  be  used  to  describe  the 
Interference  random  variables  {w( t^) } when  Impulsive  noise  (from  atmospheric 
noise  or  man-made  noise)  or  Intelligent-adversary  Interference  are  present  or 
absent.  (Note,  these  are  two  separate  problems.)  The  nominal  density  will 
describe  an  arbitrary  density  for  the  Interference  random  variables  without 
either  the  Impulsive  noise  or  Intelligent-adversary  Interference  component. 

The  contaminated  density,  occurring  with  probability  e,  will  describe  the 

* We  recognize  the  Inappropriateness  of  the  known-signal  detection  problem  model 
for  the  complete  HF  problem  ([18],  [21],  [22],  [24],  [28]).  However,  by 
treating  the  easier  robust  detection  problem  for  known  signals  first,  we  gain 
valuable  Insight  as  well  as  possible  solutions  to  the  robust  detection  problem 
for  unknown  signals. 
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interference  random  variables  containing  the  impulsive  noise  or  the  intelligent- 
adversary  interference  component.  The  interference  random  variables  will  be 
assumed  independent,  but  not  identically  distributed. 

For  this  formulation,  Huber's  results  extended  to  time-varying 
signals  ([16])  and  non-identically-distributed  random  variables 
can  be  applied  directly.  The  minimax  detector  obtained  in  this  first  step  and 
based  on  a probability  ratio  test  will  be  dependent  on  the  arbitrary 
nominal  specification.  The  next  step  will  be  to  optimize  this  minimax 
detector  over  the  class  or  appropriate  subclasses  of  possible  nominal  densities. 
The  resulting  minimax  detector  obtained  in  the  second  step  will  have  the 
robust  properties  desired  for  the  HF  communication  problem.  Steps  1 and  2 
can  be  abstractly  stated  as: 

Step  1:  Find  the  detector  and  least -favorable  mixture 
interference-density  f£  = (l-e)f0  + efj,  for  a given  e and  nominal  f0,  satisfying 


4-- 


cujff.l.  ■ *■„  f"aexF  c[,.ft] 


re  ’ e— ' c 

where  C[<f>,fe]  Is  the  probability  of  error  using  detector  <f>  and  mixture 

interference  density  fe,  0 is  the  class  of  detectors  based  on  probability 

ratio  tests,  and  Ff  = { f e | f e = (l-e)f0  + cfj;  fj  any  density}. 

Step  2:  Find  the  detector  <j>°°  and  least-favorable  mixture 

Interference-density  f®°,  for  a given  e and  class  of  possible  nominal 

densities  F0,  satisfying 

nn  min  max  _ _ 

«♦?.»?]  ' ♦°«D? 

where  C[$°,f?]1s  the  probability  of  error  using  detector  <p°  (minimax  for 

L t 

given  e and  some  f#«F0)  and  mixture  Interference  density  f°  (least  favorable 
for  given  e and  some  f#cF# ),  D°  * Is  a minimax  detector  with  respect  to 

Fe  for  given  e and  some  f^*F#},  and  F°  « (f°|f°  ■ (l-e)f.  + Fi;Vf,«F#). 
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After  evaluation  of  this  robust  detector,  additional  methods  for 
improving  performance  may  become  apparent.  For  example,  including 
an  adaptivity  capability  for  following  the  time  occurrence  of 
the  contaminated  density  would  be  advantageous  presumably  since  the 
detector  design  depends  explicitly  on  the  occurence  rate  t.  As  an 
alternative,  keeping  the  detector  non-adaptive  and  seeking  a detector 
that  would  be  robust  also  on  a range  of  e may  be  desirable.  This 


alternative  method  could  constitute  step  three,  abstractly  stated  as: 

Step  3:  Find  the  detector  d>°o * an<*  least-favorable  mixture 

interference-density  f°2,  for  a given  class  of  possible  nominal  densities 

e 

FQ  and  given  range  0 < e < e,^  satisfying 

min  max 


cr-C-O  - 


1.00^  00 


max 


foor  co 


max 


oo  oo 

CUf  .fe  ] 

*•1  2 


where  C[<f>°°,f£°]  is  the  probability  of  error  using  detector  <|>°0  (robust  for 

1 2 Ei 

some  FEl  with  0 < e,  <_  ^x)  and  mixture  Interference  density  f00  (least 

0 

favorable  for  some  D°5  with  0 < c2  < e^),  DMX  = {<t>°°|0  i ei  i Emjly}.  and 


- i - max 


- <f2!lo<  e,  i W- 


This  first  approach  should  produce  promising  results  to  apply 
to  the  HF  communication  problem.  A major  weakness  Is  assuming 
Independence  for  the  Interference  random  variables.  Consequently,  the 
performance  of  this  detector  for  dependent  Interference  random  variables 
should  be  evaluated,  probably  by  appropriate  tests  and  simulations; 
analysis  may  be  possible,  however,  for  some  forms  of  independency,  e.g., 
state  variable  models  or  Markovian  requirements. 
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Approach  2 

For  the  second  approach,  an  attempt  should  be  made  to  solve  directly  the 
formulated  problem  using  generalized  Lagrange-multlpller  theory  on  normed 
linear  spaces  ([15]).  By  modeling  the  possible,  Interference  random-sequence 
densities  as  members  of  a subset  of  a space  of  multivariate  functions,  the 
property  of  dependency  among  the  random  sequences  components  can  be  treated 
In  a convenient  manner:  you  are  not  required  to  specify  a dependency  relation- 
ship for  ,the  subset  of  multivariate  functions.  Characterizing  this  subset 
will,  however,  require  the  statistical  descriptions  of  the  Interference 
components  w^(tk),  "ou^k^’  and  d^scusse<^  *n  the  HF  comnunl cation 

problem  formulation.  This  approach  can  be  expressed  more  concisely  as: 

Step  1:  Find  the  detector  <fr°  and  least-favorable  multivariate 
density  f°  satisfying 

C[*°.  f°)  - jjj  JJJ  CU,  f] 

where  CU,  f]  is  the  probability  of  error  using  detector  0 and  Interference 
multivariate  density  f,  D Is  the  class  of  detectors  based  on  probability  ratio 
tests,  and  F Is  the  subset  of  admissible  multivariate  functions  for  which 
f * f(w)  F Implies  f Is  a probability  density  for  the  Interference  random 
sequence  » = V »QU  + "IA  * Mti). w(tN)). 

Characterizing  the  subset  of  admissible  densities  F will  probably 
be  the  most  difficult  aspect  of  this  approach  since  f(w)«F  results  from 
convolving  the  Individual  probability  densities  for  the  background  noise  Wg, 

Impulsive  noise  Wj,  other-user  Wgy,  and  Intelligent-adversary  v»IA  Interference 
components.  The  noise  model  of  Spaulding  and  Middleton  [26]  may  be  applicable 
here.  In  choosing  the  normed  linear  space,/',  for  which  Fc/'',  the  most  likely  choice 

II 

Is  Li[ft  ],  the  normed  linear  space  of  absolutely  Integrable  functions  on  the 
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domain  n ([15]).  In  this  space,  all  non-negative  functions  that  have  a norm 
of  one  represent  probability  densities. 

Another  difficult  aspect  of  this  approach  may  be  in  specifying 
exactly  the  least-favorable  density  after  applying  the  Lagrange  multiplier 
theory.  The  theory  provides  necessary  — but  generally  not  sufficient  — 
specifications  defining  the  class  of  least-favorable  density  candidates;  also  it 
provides  the  validation  method  for  identifying  the  least-favorable  density.  In 
many  instances  the  class  of  density  candidates  are  infinite  and  the  least-favorable 
density  may  be  quite  time-consuming  to  find  without  some  assistance  from 
helpful  hints  and  luck.  One  technique  for  gaining  this  investigative  assistance 
(and  other  insights)  for  multistage  problems  is  to  solve  a simpler  corresponding 
single-stage  problem.  Another  technique  is  to  solve  a more  general  problem  for 
which  the  desired  problem  is  a special  or  restricted  case.  These  techniques  will 
be  used  in  this  second  approach. 

This  second  approach  is  more  ambitious  than  the  first  approach  and 
should  yield  more  realistic  results.  The  performance  of  the  derived  robust 
detectors  in  the  actual  HF  environment  depends  on  how  closely  the  description 
of  the  subset  of  admissible  densities  F models  the  actual  environment.  If  the 
direct  approach  appears  too  formidable  in  this  regard,  this  closeness  of  the 
model  to  the  actual  environment  may  be  accomplished  in  successive  steps.  The 
first  step  (Step  la)  could  be  equivalent  to  Step  1 above  with  ¥'  defined  as  a 
larger  subset  of  $ : for  example,  the  subset  containing  all  probability  densities. 

A valid  corresponding  second  step  (Step  lb)  could  be  a sensitivity  investigation 
of  the  Step  la  robust  detector  subjected  to  deviations  from  the  corresponding 
least-favorable  density.  These  deviations  would  be  typical  densities  from  the 
class  F defined  in  Step  1 above.  Hopefully,  some  least-favorable  deviation 
could  be  found  in  this  step.  Another  valid  corresponding  second  step  (Step  lb') 
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would  be  equivalent  to  Step  la  with  F'  replaced  with  a proper  subset  F"c:F', 
l where  F<=f".  This  series  of  nested  problems  (Step  la.  Step  lb'  , . . .)  can 

continue  In  principle  until  the  problem  Step  1 Is  reached  or  approximated 
closely  enough  for  practical  purposes.  An  advantage  of  this  nested  problem 
approach,  hopefully  for  little  effort  in  characterizing  F',  F",  ....  is  the 
considerable  Insight  likely  gained  as  to  how  the  corresponding  robust 
detectors  <|>0',<i>0,/,  ....  approach  the  robust  detector  <p°  of  Step  1. 
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